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FOREWORD 


This report covers the final work completed on the research project 
"Effects of Angle of Attack on the Coupled Radiative and Convective Heat 
Transfer About Blunt Planetary Entry Bodies." The work was supported by 
the NASA/Langley Research Center (Aero thermodynamics Branch of tha Space 
Systems Division) through research grant NSG-1464. The grant was Initially 
monitored by Dr. Randolph A. Graves, Jr. of the Space Systems Division. 

In recent months, the grant has been monitored by Mr. Jim J. Jones of the 
Spac«. System Division. 
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LAMINAR AND TURBULENT FLOW SOLUTIONS 
WITH RADIATION AND ABLATION INJECTION 
FOR JOVIAN ENTRY 


By 

AjAy Ktaar 1 and S. N. Tiwari 2 
SUtMARY 

Laminar aid turbulent flow- field solutions with coupled carbon-phenolic 
■ass injection are presented for the forebody of a probe entering a nominal 
Jupiter atmosphere. Solutions are obtained for a 35-degree hyperboloid and 
for a 45-degree spherically blunted cone using a tine-dependent , finite- 
difference aethod. The radiative heating rates for the coupled laminar flaw 
are significantly reduced as compared to the corresponding no-blowing case; 
however, for the coupled turbulent flow, it is found that the surface radiative 
heating rates are substantially increased and often exceed the corresponding 
no-blowing values. Turbulence is found to have no effect on the surface 
radiative heating rates for the no-blowing solutions. The present results 
are compered with the other available solutions, and soma additional solutions 
e.*e presented. 


INTRODUCTION 

The entry probe for the NASA Galileo mission to the planet Jupiter will 
encounter very severe aerodynamic heating due to exceedingly high entry 
velocities. A massive heat shield will be required on the probe's surface 
to protect the scientific instruments from this aerodynamic heating. To 
reduce the uncertainties in the design of the heat shield. It is necessary 
to accurately predict the entry aerothermal environment. Since it Is 
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extremely difficult to duplicate the entry enviro n ment in an experimental 
facility, the design of the heat shield for the Galileo probe relies primarily 
on analytical predictions. 

In order to predict the flow-field accurately, one must modal all the 
physical phenomena occurring in the flow. For the Galileo proble, the gases 
surrounding the probe forebody will have temperatures In the range of 10,000 
to 20,000 K, and the surface heating will be primarily due to radiative 
heating. This will result in massive ablation of the probe's hear shield. 
Moreover, the flow is expected to be mostly turbulent. The analytical 
technique should include all the above-mentioned phenomena in order to 
accurately predict the surface heating rates sod other flow characteristics. 

Reference 1 presented laminar and turbulent flow-field solutions for 
the analytic-shaped bodies (hyperboloids) using a nominal entry trajectory 
for Jupiter. The viscous shock-layer analysis included coupled ablation 
injection, detailed spectral calculation of radiative heating, and 
equilibrium chemistry. In reference 2, this analysis was further v«ed to 
obtain some coupled turbulent flow-field solutions for the forebody of a 
45-degree spherically blunted cone (Gallleo-probe configuration) only at 
certain times of the nominal entry trajectory. No solutions could be obtained 
for the 1 — coupled flow field. Reference 3 has presented coupled 
turbulent solutions for the Galileo probe using an approximate nethod which 
assumes a shock shspe and then predicts the body shape. There still is a 
great need for an analytical technique which can accurately predict the 
coupled laminar and turbulent flow fields for the trajectory. The purpose 
of the present study was to develop such a technique. 

In reference 4, a time-asymptotic, finite-difference method was developed 
to predict the forebody flow field past blunt axlsyMetric bodies, such as 
spherically blunted cones, hyperboloids, etc., at zero and small angles of 
attack. The method uses time-dependent, viscous, shock- layer equations to 
describe the flow field. The analysis was modified in reference 5 for the 
laminar and turbulent flow of a radiating and reacting gas under chemical 
equilibrium. It was found In reference 6 that the method works well for both 
the laminar and turbulent flows even In the presence of massive surface 
blowing. 
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The present investigation um the above method for obtaining the 
laminar and turbulent flow-field solution* for a 35-degree hyperboloid and 
a 45-degree spherically blunted cone (Galilee probe shape). The analysis 
includes coupled ablation injection, detailed spectral calculation of 
radiative heating, and equilibria chemistry. The radiative heat transfer 
and equilibrium chemistry calculations are made using the sane techniques 
used in reference 1. The results of the present analysis are compared with 
those of references 2 and 3, and additional solutions are presented for 
the conditions for which reference 2 presents no solution. 

LIST OF SYMBOLS 



ablator mass fraction 

mass fraction of species i, p^/p 

mass fraction of element i 

total enthalpy, H/V? 

V 

enthalpy of mixture, e C h 

i - 1 11 

enthalpy of undecomposed ablation material 

enthalpy of species i, h^/V^ 

thermal conductivity of mixture, k^y^V^/T^) 

Lewis number 

turbulent Lewis number 

molecular weight of mixture 

molecular weight of species i 

molecular weight of element l 

mass injection rate, 

number of species 

coordinate normal to the body, n/R 

n 

number density of atomic hydrogen (particles/cm 3 ) 

— — o 

pressure, p pJT 
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net radiative heat flux In n-dlreetlon, q /p V’ 
component of radiative flux toward* the shock 
component of radiative flux towards the wall 
convective heat flux to the tjall 
universal gas constant (J/aol*°K) 
nose radius (■) 

free-streaa Reynolds nimber , p V R /u 

• °* o • 

radium measured from axis of symswtry to point on body surface, r/R 

D 

coordinate measured along the body, s/R 

H 

Schmidt number 
temperature, T/T 

Oft 

free-stream temperature (°K) 

nondimen8ional time, tV /R 

• n 

tangential velocity, u/V^ 
tangential velocity at shock, u^/V^ 
normal velocity, v/V m 
free-stream velocity (m/sec) 
r + n cos 0 

coefficients for sublimation temperature [eq. (5)] 

shock standoff distance, 5 _/R_ 

number of atoms of the Zth element in species 1 

emlsslvity 

normalized eddy viscosity, u^/u 
body angle measured from the body axis 
local body curvature, icR n 
1 + nic 
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u viscosity, Zh m 

Ry addy viscosity, 

u m free-atream viicoilty (N ••c/m 2 ) 

p density, ~c/o m 

0 m fr«t>ttrMi density (kg/u 3 ) 

a Frandtl mabtr 

Oj turbulent Prandtl number 

a Stefan Boltzmann constant 
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Superscript: 

( ) dimensional quantities 
Subscript 

1 1th species 

t l th element 

s shock values 

sub subllnatlon 

v vail values 

• free-stream values 

— values for the solid ablation material at the surface 

o stagnation point 



ANALYSIS 

Governing Equations 

The equations of motion for the laminar and turbulent flows of reacting 
gas mixtures In chemical equilibrium are given in reference 5. These viscous 
ahock-layer-typs equations are modified for the present analysis to include 




5 








ch« elemental continuity equation which la required here due to the surface 
ablation injection. The governing equations in the body-orl*oted coordinate 
systea (fig. 1) for the flow at saro angle of attack are expressed as 


where 
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The equation of state is given by 


P« 



( 3 ) 


Two independent variable transformations are applied to the governing 
equations. The first transformation maps the computational domain into a 
rectangular region in which both the shock and the body '< -e made boundary 
mesh lines. The second transformation further maps the computational region 
into another plane to allow higher resolution near the body surface without 
any significant increase in the number of mesh points in the normal direction. 
The details of these transformations are given in reference 4. 

Boundary Conditions 

No-slip boundary conditions are used at the surface. The wall temperature 
and mass injection rate are either specified or calculated. The boundary 
conditions at the shock are calculated by using the shock relations. 

For the calculated mass injection conditions, the ablation process is 
assumed to be quasi-steady and the wall temperature is the sublimation 
temperature of the ablator surface. The coupled mass injection rate is then 
given by 
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The sublimation temperature for the carbon-phenolic ablator is given by 
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ORIGINAL jKdui'j,b\ 

OF POOR QUALIFY 



where is the vail pressure In atmospheres and is the ablator mass 

fraction at the wall. The values of 6 are given in table 1. These 

J 

are applicable for a free-stream gas composition of 89 percent H 2 * 11 

percent He by volume, and for ablator mass fractions of 0.4 to 1.0. 

For ablation injection, the elemental concentrations at the wall are 
governed by convection and diffusion as given by the equation 



£ 

The net r-diative flux, q , can be represented as 


R 

q * 



R~ 
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(7) 


At the surface 



a T 1 * 
w 


( 8 ) 


The surface is treated as a gray surface with an assumed emissivity e of 
0.8 and a reflectivity of 0.0. The heat transferred to the wall du« to 
conduction and diffusion is 


C 

-«w 



N 

r 

i - 1 



(9) 


Radiative Transport 
R 

The radiative flux, q , is calculated with the radiative transport 
code RAD (refs. 7,8) which has been incorporated into the present computer code. 

The RAD accounts for the effects of nongray self absorption. Molecular 
band, continuum, and atomic line transitions are included. A detailed frequency 
dependence of the absorption coerfcients is used for integrating over the 
radiation frequency spectrum, and the tangent slab approximation is used for 
integrating over physical space. The chemical species considered in the 
present study for determining the radiative transport are H, H 2 , H + , e , C, C 2 , 
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Cj, C , C , CO, 0, O 2 , 0 , and 0 . The radiation properties used in this 
study for C 2 and C 3 are those reported in references 9 and 10, 
respectively. 


Thermodynamic and Transport Properties 

The equilibria!, composition is determined by a free energy minimization 
calculation as developed in reference 11 and written for the present code in 
reference 12. Thermodynamic properties for specific heat, enthalpy, and 
free energy and transport properties for viscosity and thermal conductivity 
are required for e_uh species considered. Values for the thermodynamic and 
transport properties are obtained by using polynomial curv<r| fits. Mixture 
viscosity is obtained by using the semi-empirical formula of Wilke (ref. 

13). 

Six chemical species are used to describe the hydrogen-helium gas 

mixture: e , H, H 2 , H , He, and He . With carbon-phenolic injection, 13 

+ + 
additional species are used: C, C 2 , C 3 , C , C 2 H, C 3 H, C 4 H, C 2 H 2 , 0, O 2 * 0 , 

CO, and CO 2 . The Lewis number and Prandtl number of the mixture are set 

equal to 1.1 and 0.64, respectively. 

Eddy-Viscosity Approximations 

A two-layer, eddy-viscosity model consisting of an inner law based upon 
Prandtl' s mixing-length concept and the Clauser-Klebanoff explosion (based on 
refs. 14 and 15) for the outer law is used in the present investigation. This 
model. Introduced by Cebeci (ref. 16) , assumes that the inner law is applicable 
for the flow from the wall out to the location where the eddy viscosity given 
by the inner law is equal to that of the outer law. The outer law is then 
assumed applicable for the remainder of the viscous layer. It is noted that 
the eddy viscosity degenerates to approximately zero in the inviscid portion of 
the shock layer. The degeneracy is expressed in terms of the normal intermit tency 
factor given by Klebanoff (ref. 15). Reference 2 gives a more detailed 
description of the turbulence model and various expressions for it. The only 
difference between the models used in the present calculations and that used 
in reference 2 is in the boundary- layer edge definition. The present analysis 
uses the definition given in reference 17, which is based on an index of 
diffusion, conduction, and dissipation. The turbulent Prandtl number and Lewis 
number are assumed to be 0.9 and 1.0, respectively. 


'V 
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METHOD OF SOLUTION 


A time- asymptotic, two-step, finite-difference method due to MacCormack 
(ref. 18) is used to solve the governing equations. The details of the 
method are given in reference 4. 

The calculations of equilibrium chemistry, radiative heat flux, and 
eddy viscosity require a significant amount of computing time, and it is not 
feasible to perform these calculations in each time-step. In the present 
analysis, the eddy viscosity is calculated after every 25 time-steps, the 
equilibrium chemistry after every 200 time-steps, and the radiative heat 
flux after every 1,000 or 2,000 time-steps. 

RESULTS AND DISCUSSION 

Forebody flow-field solutions and surface heating rates are presented for 
a probe entering a nominal Jupiter atmosphere. Most of the results are 
obtained for the flow past a 45-degree sphere-cone (the current Galileo 
probe configuration) , but a sample calculation is made for the flow past a 
35-degree (asymptotic half angle) hyperboloid for comparison with the 
existing results. The free-stream conditions correspond to ar, entry 
trajectory into the Jupiter atmosphere where the atmospheric gas model is 
the Orton nominal atmosphere consisting of 89 percent of H? and 11 percent 
Re by volume (see ref. 2 for further details). These free-stream conditions 
are taoulated in table 2. The nose radius of the probe is 0.3112 m. 

In this study, solutions are obtained for both laminar and turbulent 
flows with and without coupled ablation injection. In all the turbulent 
solutions, transition from laminar to turbulent flow is assumed to occur 
Instantaneously at the first mesh point away from the stagnation point, 
which is at s ■ 0.157. Both laminar and turbulent solutions are compared with 
the existing results for the hyperboloid, but only turbulent solutions are 
available for the sphere cone. 

Figures 2 and 3 show the surface heating rates with and without coupled 
ablation injection for the sphere-cone at 111.3 sec. It is seen from figure 2 
that the turbulence has no effect on the radiative heating rates for no-injection. 
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but the turbulent convective heating rates are much higher than the laminar 
convective heating rates, as is seen from figure 3. The radiative heating 
rates for the coupled laminar flow are significantly reduced all over the 
body and the convective heating rates are negligible. However, the radiative 
heating rates for the coupled turbulent flow are seen to Increase again as 
compared to the coupled laminar flow and, in fact, are even higher than the 
no- Injection radiative heating rates on the conical portion of the probe. 

This adverse effect of turbulence on the radiative heating rates was first 
reported in reference 2. The coupled turbulent convective heating rates also 
increase but still are only a small fraction of the corresponding radiative 
heating rates. Thus, it is concluded from figure 2 that the coupled turbulent 
radiative heating rates are critical for the design of the probe's heat shield. 
The enhancement of the radiative heating rates due to turbulence will be 
discussed later. 

Figure 4 shows a comparison of coupled radiative heating rates for the 
laminar and turbulent flows over the 35-degree hyperboloid with reference 2. 

It is seen that the present laminar radiative heating rates are about 5 to 
10 percent higher and the present turbulent radiative heating rates are about 
20 percent higher than those predicted in reference 2. This difference is 
partially due to the fact that, in reference 2, the mass injection rate at a 
particular body station is calculated by using the heating rate at the previous 
body station, wince the heating rates over the body are seen to decrease 
with increasing body distance, the mas* injection rate at a particular body 
station used in reference 2 is higher than it should be if the heating rate at 
that same body station is used. This higher mass injection rate results in 
a slight decrease in the radiative heating rate. Ths present analysis uses 
the heating rate at the same body station at which the mass injection rate 
is being calculated. Another difference between the present results and those 
of reference 2 is that the present analysis does not predict any overshoot in 
the turbulent radiative heating rates near the stagnation region. 

Figure 5 shows the comparison of the corresponding surface mass Injection 
rates. It is seen that although the present laminar radiative heating rates are 
slightly higher than those in reference 2, the present surface mass injection 
rates are still lower. This is directly due to the reason given earlier. The 
present turbulent mass injection rates are higher on the downstream part of 
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the body, but the difference Is not as much as it was between the radiative 
heating rates. This decrease in the difference in the mass injection rates 
is also due to the reason discussed earlier. 

Shock standoff distances are compared in figure 6. The present analysis 
predicts about 10 percent high-- shock standoff distances for the turbulent 
flow, but the shock standoff distances for the laminar flow compare well. 

All the results presented below are for the flow past a sphere cone. 

Figures 7 and 8 show the coupled laminar and turbulent radiative heating 
rates and surface mass injection rates at 111.3 sec. The results of reference 
2 are also plotted for comparison. It is seen that the present results are 
about 10 to 15 percent higher on the conical portion of the probe. Unlike 
reference 2, the present heating rates and mass injection rates gradually 
decrease on the nose portion of the probe. No comparison could be made for 
the coupled laminar flow as the method of reference 2 does not work for the 
coupled laminar flow past a sphere cone. 

Figure 9 shows the surface pressure distribution for the coupled laminar 
and turbulent flows at 111.3 sec along with the surface pressure distribution 
predicted in reference 2 for the coupled turbulent flow. It is seen from the 
present results that the surface blowing smoothes out the effect of curvature 
discontinuity at the juncture point of the probe. The surface pressures 
for the turbulent flow as predicted by the two analyses compare well in the 
nose region and on the conical portion of the probe, but the results of 
reference 2 show a sharp decrease and then Increase in the pressure around 
the juncture point. This behavior of the surface pressure in the presence 
of massive surface blowing is not physically expected. 

Figures 10 and 11 show the coupled laminar and turbulent radiative 
heating rates at 109 and 107.2 seconds. The method of reference 2 fails to 
work for these trajectory points, but coupled turbulent radiative heating 
rates are given in reference 3. The present radiative heating rates are 
about 25 percent higher than those predicted in references 3. Figure 12 shows 
the corresponding mass injection rates. The nondimenslonaJ mass injection 
rate m is maximum for 107.2 sec out of the three trajectory points 

e 

considered in the present analysis. For other trajectory points, m is also 
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expected to be lower than for 107.2 aec. Since the present method worked well 
for 107.2 sec. It should be able to provide solutions at any point on the 
trajectory given In table 2 . 

Figures 13 through 16 show the effect of turbulence on the shock layer. 

It Is seen from these figures that the turbulence significantly changes the 
structure of the shock layer. The profiles of reference 2 for the turbulent 
flow are also plotted for comparison. In general, the present flow profiles 
show distributions similar to those in reference 2. Tangential velocity and 
temperature profiles are shown in figures 13 and 14. These profiles show 
the same behavior as observed in reference 6 in the presence of massive 
surface blowing. The gradients near the surface significantly increase for 
the turbulent flow and t.ie turbulence brings the high-temperature gases closer 
to the surface. Figure 15 shows that, for the turbulent flow, the ablation 
layer (the portion of the shock layer in which the ablation products are 
dominantly present) is much thicker, but the concentration of the ablation 
products is decreased due to enhanced diffusion. Figure 16 shows the number 
density distributions of Cj and C 3 for laminar and turbulent flows. It 
is seen that the number densities of C 2 and C 3 for the turbulent flow are 
significantly reduced as compared to the corresponding laminar values. 

The high temperatures near the surface for the turbulent flow dissociate these 
molecules, thus reducing their concentrations. It is shown in reference 10 
that Cj and C 3 are the main radiation-absorbing molecules for Jovian entry 
conditions. Thus, the turbulence has an adverse effect on the effectiveness 
of the ablation species to absorb radiation. 

CONCLUSIONS 

A time-dependent, finite-difference method is used to develop a code for 
solving Che coupled laminar and turbulent flows over the forebody of a probe 
entering a nominal Jupiter atmosphere. Solutions are obtained for both a 
35-degree hyperboloid and a 45-degree sphere cone. Detailed comparisons are 
made with the existing results. It is found that the present code works well 
for both the coupled laminar and turbulent flows at all the trajectory points. 
The code used in reference 2 works well for both the coupled laminar and 
turbulent flows over a hyperboloid, but it works only for the coupled turbulent 
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flows ovar a sphere cone. Even for the coupled turbulent flows, it does 
not work at several trajectory points where the nondimens lonal Injection 
rates are very high, such as at 109 or 107.2 sec. The present coda uses 
the heating rate at the same point for which the injection rate is being 
calculated, whereas in reference 2 the heating rate at the previous mesh 
point is used. 

The turbulence is found to alter the structure of the shock layer in 
such a way that the effectiveness of the ablation layer is significantly 
reduced. The present analysis supports the finding of reference 2, in which 
the coupled radiative heating rates exceed even the no-injection values on 
the conical portion of the body. The present coupled, turbulent, radiative 
heating rates have the maximum value at the stagnation point, whereas the 
results of references 2 and 3 show the maximum to be slightly downstream of 
the stagnation point. The surface pressures predicted by the present analysis 
appear to be more physically plausible in the presence of massive surface 
blowing than those in reference 2, where a very sharp dip occurs around the 
sphere-cone juncture point. Detailed flow-field profiles are similar to 
those in reference 2, but there are differences in the magnitude. 



Table 1. Sublimation temperature coefficients 


j 


Coefficients 



1 

2 

3 

4 

5 

6 i.J 

5552 

-20184 

53058 

-57933 

23243 

8 ’.J 

1798 

-12049 

30145 

-32045 

12457 

*3.3 

322 

- 2208 

5270 

- 5450 

2092 


Table 2. Free-stream conditions. 


Time from Altitude * ' « 

entry point (km) (km/sec) (kg/m 3 ) 
(sec) 


99.9 

195.56 

103.9 

169.57 

107.2 

149.10 


138. 60 

110.2 

131.90 

111.3 

126.05 

112.2 

121.49 

113.5 

115.31 

114.2 

112.20 

115.3 

107.63 

116.4 

103.45 

117,4 

100.00 


48.07 

2. 65*10’ 5 

46.96 

7. 19*10’ s 

44.83 

1.64*10’ 4 

42.88 

2. 54*10’ 4 

41.16 

J.38*10’ 4 

39.29 

4.36*10’ 4 

37.52 

5. 34*10" 4 

34.67 

7.02*10’ 4 

33.01 

8.05*10’ 4 

30.31 

9. 89*10’ 4 

27.54 

1.20*10’ 3 

25.07 

1.41*10’ 3 
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Figure 2. Effect of coupled carbon-phenolic injection on radiative 
heating rates for a 45-degree sphere cone at 111*3 sec* 


19 




300 


LAMINAR TURBULENT 



Figure 4. Comparison of coupled radiative heating rates for a 
35-degree hyperboloid at 111.3 sec. 
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Figura 7. Comparison of coupled radiative heating rates for a 
45-degree sphere cone at 111.3 sec. 
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Figure 9. Comparison of surface pressure distribution for a 45-degree 
sphere cone at 111.3 sec. 






Figure 10. Comparison of coupled radiative heating rates for a 45-degree 
sphere cone at 109 sec. 
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Figure 13. Comparison of turbulent tangential velocity profiles at 
s - 1.8 for a 45-dep,ree sphere cone at 111.3 sec. 
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Figure 16. Effect of turbulence on concentrations of major radiation 

absorbers at s m 1.8 for a 45-degree sphere cone at 111.3 sec. 
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